Geometric phase in open systems 
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We calculate the geometric phase associated to the evolution of a system subjected to decoherence 
through a quantum-jump approach. The method is general and can be applied to many different 
physical systems. As examples, two main source of decoherence are considered: dephasing and 
spontaneous decay. We show that the geometric phase is completely insensitive to the former, i.e. 
it is independent of the number of jumps determined by the dephasing operator. 
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The generalization of the geometric phase 0, 0, as ~ 
sociated to the evolution of a system in a mixed state sce- 
nario is still an open problem. There have been various 
proposals tackling the problem from different perspec- 
tives like, for example, via state purification p|, through 
an interferometric procedure ||, or for dissipative sys- 
tems 0. In most of the cases these definitions do not 
agree on account of different constraints imposed, namely 
different generalizations of the parallel transport condi- 
tion. 

In this paper, we approach the problem of mixed states 
geometric phases in open quantum systems through the 
quantum jumps method 0, 0] (for alternative methods 



see UJ, [10|, lUl). We show that, in general, it is possi- 
ble to define and calculate a geometrical phase for an 
open system, which can even be robust against decoher- 
ence, as, for example, it is the case for purely diffusive 
reservoirs. As an illustration of our method, we calculate 
Berry's phases for a spin 1/2 in contact with different 
reservoirs. 

The quantum jumps model proves to be particularly 
suitable for the case of geometric phases because in each 
particular trajectory the quantum state of the system re- 
mains pure (if initially pure). The non unitary evolution 
of the system and mixed states are recovered when aver- 
aging over all possible trajectories. Therefore, we avoid 
the problem of finding the proper parallel transport con- 
dition, which is very well defined for unitary evolutions. 

Let us take a system evolving according to the follow- 
ing general master equation (H = 1): 

1 1 ™ 

p = -[ H ^p\-^Y,{ T l T kp + p^k -2r feP r+}, (1) 

fe=i 

For a small time interval At, we can describe the time 
evolution of the density matrix by 



p(t + At)*Y,W kP (t)Wl, 



(2) 



fe=0 



where W = 1 - iHAt and W k = VAtT k (k € 
{1 . . . n}) are called the "no-jump" and jump operators 



respectively, 
by: 



H is a non-Hcrmitian Hamiltonian, given 



L k=l 



(3) 



Note that the operators W k fulfill the completeness rela- 

In this description, the dynamics of the system is ap- 
proximated by dividing the total evolution time T into a 
sequence of discrete intervals At = ■£>. According to 
Eq. (J2J), the state of the system, after any time step 
t m = mAt, evolves into p(t m+ i) = Wkp(t m )Wl (up to 
first order in At), with probability p k = TrWkp{t m )Wl. 
For example, an initial (t = 0) pure state ipo would 
evolve, after the first time interval, into the (not nor- 
malized) state \ipo) — ► IV'i) = Wk\ipo) W1 th probability 

Pfc(ti) = (Vto|W+Wfchto>. 

The time evolution of the system is, then, calculated 
for a set of possible trajectories containing, each one of 
them, different numbers of jumps, occurring at different 
times, i.e. each trajectory is defined as a chain of states 
obtained by the action of a sequence of operators Wk on 
the initial state. For example, for an initial pure state 
IV'o); the (non-normalized) state of the system, after the 
m-th step, along the i-th trajectory, is given by: 



U.w m \iM 



(4) 



i=i 



where stands for the 1-th element of a sequence of 
indexes with values belonging from . . . n. Each trajec- 
tory, then, is represented by a discrete sequence of pure 
states {if)a, tp®, ■ ■ ■ , tp^}- The dynamics given by the 
master equation is recovered by summing incoherently 
all the states associated to each trajectory, and taking 
the continuous limit At — > 0. 

The fact that a pure state remains pure in each trajec- 
tory, in the quantum jumps method, is very useful in our 
case, since it is known that, given a chain of pure states 
. . .\iPn}}, the geometric phase associated to them 
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is given by the Pantcharatnam formula lz 



la 



;{ (V'iI^)^!^) ■ • ■ (^n-i\M(M^i)} (5) 



Therefore, we are able to associate a meaningful geomet- 
rical phase to each trajectory "i" described by the sys- 
tem, as the continuous limit of Eq.© for the sequence 

As an example, let us consider the "no-jump" trajec- 
tory for a completely general master equation. The evolu- 
tion of a quantum state along this trajectory is obtained 
by the repeated action of the operator Wq. At the time 
t = mAt, the quantum state will be approximately given 
by: 



10 - (w rm = 



l-i—H 

N 



IVto). (6) 



which in the continuous limit N — > 00 yields to a dy- 
namics governed by the complex effective Hamiltonian 
H: 

ij t \^°(t))=H\^(t)) |V>») = h/> ) (7) 

Thus, the evolution corresponding to this trajectory is 
given by a smooth chain of (non normalized) states \ip(t)), 
in which case 7jv converges to: 



7 = —Im 



T W)\&Mt)) 
W)W)) 



dt-arg{(^(T)|V(0))}. (8) 



Substituting Eq. Q into Eq.(jHJl, we obtain the geometric 
phase for a no-jump trajectory, which is given by: 



(Tp°(tW(t)) 



di-arg{(^(T)|V>°(0))} (9) 



This is the geometric phase associated to a non-unitary 
evolution of a system |12l Il3j]. when there are no jumps. 
The first term is clearly the opposite of the dynamical 
phase associated to the non-unitary evolution, as it is 
given by the average of the Hamiltonian (up to a minus 
sign) along the path traversed by the system. The second 
term is the total phase difference between the final and 
the initial state, according to Pancharatnam's definition 
of distant parallelism Thus the geometric phase is 
obtained as the difference between total and dynamical 
phase associated to a given evolution of pure states 

Note that, in the special case in which ^l^i ^ 1 

(which is a unital evolution) , the geometric phase associ- 
ated with the no-jump trajectory is the same as the one 
acquired by an isolated system evolving under the same 
Hamiltonian H . This becomes clear when one notes that, 
in this case, Wq = (1 — a)l + iHAt and the evolution 
of state \^(t)) is the same as its isolated counterpart up 
to a global normalization factor e~ at . In other words, 
for this particular source of decoherence, if the reservoir 



is permanently measured and no jump is detected, there 
is no gain of information on the system, which simply 
projects it back into its unitary evolution. 

Note, also, that following the idea of [l4|, it is possible 
to represent the geometric phase (j5J) as the integral of 
the Berry connection form: 



doj = Im 



(m\m) 



(10) 



along a closed path. This path is formed by the trajec- 
tory i/j(t) followed by the states along the Hilbert space 
during the dynamical evolution and the shortest geodesic 
connecting final and initial states ip(T) and ^(0). Thus 
the second term of equation can be regarded as the 
path integral of the Berry connection along this geodesic. 

Suppose, now, that there is only one jump in the tra- 
jectory at an arbitrary time t\, which occurs in a time 
much shorter than any other characteristic time of the 
system. Then, we can separate the evolution in two parts 
(before and after the jump) and, the continuous limit of 
equation JSJ) leads to the following expression: 
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dt- a vg{(^ (M|r,# (0))} + 

dt-arg{<v"(T)|^'(0))} (11) 



where Wj is the operator associated to the occurred 
jump, and ip (t) and ip (t) are the states evolving under 
the effective Hamiltonian H, before and after the jump 
respectively. They are given by the equation JJJ) with 
initial conditions ip (0) = Vo and "0 = Wjtp (ti), 
respectively. 

The first and third term represents the dynamical 
phase given by the effective evolution Q, before and af- 
ter the jump occurs. The last term is the phase difference 
between initial and final state of the total evolution. The 
second term is a phase associated to the occurrence of 
a jump at time t\. Analogously to the total phase as- 
sociated to final and initial state, this term represents 
the phase difference between the states after and before 
the jump, and geometrically, it can be regarded as the 
path integral of the Berry connection along the shortest 
geodesic joining them. 

This result can be easily generalized to any trajectory, 
allowing for a more complicate sequence of jumps and 
no-jump evolutions. The geometric phase is then repre- 
sented as the sum of terms of the form arg (tp(ti)\Tj\ip(ti)) 
regarded as the phase associated to the jump Fj occur- 
ring at the instant ti, and terms of the form (JSJ) for the 
no jump evolutions. And clearly all these phases can 
be regarded as the integrals of Berry connection along a 
complex path composed of geodesies joining initial and 
final state of the jumps, and the paths traversed by the 
state during the evolution under H. 
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Let us apply this general quantum jumps procedure 
to a well known physical system. First, let us consider 
the simplest example of decoherence: a two levels sys- 
tem evolving under the free Hamiltonian H = %<J Z and 
subjected to dephasing, which can be described by the 
Master equation JIJ with T = \a z , where A is the the co- 
efficient giving the probability per unit time of a "phase- 
jump" . 

Since this is a decoherence model for which T^T = 
a z oc 1, which is a simple instance of a unital evolu- 
tion, according to the previous considerations, the geo- 
metric phase associated to the no-jump case is given by 
the standard geometric phase associated to the unitary 
evolution of a spin 1/2 linearly coupled to a constant 
magnetic field. For instance, after a time t = 2tt/u}, 
7o = 7r(l — (ipoW z\tpo)) = t(1 — cos6>), where tp is the 
initial state and 9 is its azimuthal angle in the Bloch 
sphere representation. 

Although the no-jump case may seem trivial, this sys- 
tem has a much more remarkable property: the geometric 
phase is actually robust against dephasing, in this sim- 
ple, but very useful example. In fact, we show below that 
the final geometric phase is unaffected by any number of 
jumps for any particular trajectory. To show that, Let 
us consider first the case of a single jump, in which the 
phase is given by: 

7i = - J -^(tpoWzmidt - arg{(V>oM"0o)} 

- J -^{ipoWz\ipo)dt 

- argj^ole^^^V.e^^lVo)} 
= tt(1 - (lp \<T z \lp )) = tt(1 - COS0), 

where the fact that H and T commute has been used. 
This result is easily generalized to any number k of jumps: 

Ik = - J -(ip \cr z \ijjo)dt- arg{(^ok z |-0o> fc } - 

- arg^ole^^lVo)} = tt(1 - costf), 

Thus, no matter how many jumps occur in the chosen tra- 
jectory, we can associate the same geometric evolution to 
the system. There is a simple geometrical explanation for 
this effect. Dephasing is a special source for decoherence 
because it does not change the projection of the spin vec- 
tor on the direction of the magnetic field, i.e. it does not 
change the relative angle 6 between the directions of the 
magnetic field and the spin. After each jump, the spin 
is still precessing around the magnetic field alongside the 
same curve. As a result, the total area covered by its 
trajectory remains the same, and so does the geometric 
phase acquired by the spin state, which is proportional 
to this area. Therefore, in the end, the geometric phase 
acquired by the spin state will be the same, no mat- 
ter how diffused its total phase may be. That does not 



mean that dephasing will not affect the measurement of 
this phase. Indeed, it will lower the visibility of any in- 
terference measurement made on the spin, because the 
visibility of the state is lowered when its mixedness is 
increased (we will address this in more details in a sep- 
arate publication). However, as the calculations above 
show, the reduced visibility will be caused by a random- 
ization of the dynamical phase, and not the geometrical 
one, which proves to be much more robust in this case. 

A more realistic example includes spontaneous decay 
as a source of decoherence for the spin 1/2 system. In 
this case, it is only worth analyzing the no-jump case, 
since any jump causes immediate and complete loss of 
phase information of the quantum state. Spontaneous 
decay T — a<7- is a decoherence source that cannot be 
associated to a unital map (cr_|_<7_ ^ 1) and, therefore, 
the phase will be affected even if no jump is detected. 
However, as we show in figure 2, the no-jump trajectory 
is a smooth spiral converging to the lower state, which 
still allows us to calculate the phase using Eq. JHJ. We 
obtained 7 = 7r + In ((-0 O |e — 47r ^ CTi | Vo)) j which in the 
limit uj a leads to 

7 m 7r(l - cos6») + (47r) 2 - sin 2 + o ( — ) 

Again, this result has a very simple geometrical explana- 
tion: as we observe the reservoir and detect no jump, the 
probability that the system is in the lower state smoothly 
increases, changing 6 and, therefore, the element of area 
covered by the spin trajectory in each infinitesimal time 
interval, as shown in figure 2. 

Another simple case that can be analyzed is the spin 
flip alongside an arbitrary direction Y = a 7 \ . In this case, 
the no jump situation is again trivial and similar to the 
dephasing reservoir, since er? = 1. When one or more 
jumps occur, we can use Ea. (|ll|) (or its generalization to 
many jumps) to easily calculate the final phase, which 
will be a sum of the partial areas covered in each trajec- 
tory with plus or minus sign depending on the respective 
coupling energy of the spin with the magnetic field. Our 
treatment is, of course, applicable even when the master 
equation contains many different sources of errors act- 
ing simultaneously on the system, since we can use the 
generalized form of equation (|ll|l to calculate the phase. 

In conclusion, in this paper, we present a method to 
calculate geometric phases in open systems. Our method 
is general and can be applied as long as the system dy- 
namics is described by a master equation in the form 
of Eq. , which is the most general c omp letely positive 
trace preserving continuous evolution |15| . By using the 
quantum jumps approach we avoid the problem of defin- 
ing Berry's phases for mixed states: in each trajectory, 
the quantum state of the system remais pure and the 
phase can be calculated through usual procedures. In 
particular, we show that it is always possible to calcu- 
late this phase, either for the no-jump trajectories or 
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FIG. 1: Evolution of the state along a "one-jump" trajectory 
on the Bloch sphere, in the case of phase diffusion decoherence 
(r oc <7z). From time to to time ti the state evolves under the 
no-jump hamiltonian H alongside the parallel of the sphere. 
At time ti a jump occurs, nipping (instantaneously) the Bloch 
vector about the z to the point t[, and the no-jump evolution 
starts again. At time ti = to + 2n/u the geometric phase 
7 = 7r(l — cos 6*) is recovered. The geometric phase is half the 
area enclosed in the path spanned by the Bloch vector. This 
is given by two contribution, the surface ti — to — t[ and the 
surface to — ti — ti. 



for the ones in which one or more jumps occur. We 
also show that, for special unital decoherence sources, 
the phase remains unaffected for the no jumps trajecto- 
ries. As a direct application of our method, we calculate 
the geometric phases of spin 1/2 systems coupled to dif- 
ferent reservoirs. We show that those phases are totaly 
robust against phase diffusion, in which case the lower 
visibility observed due to the non-unitary evolution may 
be attributed solely to a randomization of the dynami- 
cal phase. This property may be interesting for possible 
applications, specially in quantum computing, since de- 
phasing may be difficult to monitor and correct, in gen- 
eral. Therefore, it is interesting noticing that geometric 
phases are robust against this decoherence source. We 
also present a nice geometrical explanation to this effect, 
as well as to the effect on the geometric phase when spon- 
taneous emission is present, but no jump is detected. We 
also briefly comment on other typical decoherence effects 
on the system, like arbitrary spin flips. The method pre- 
sented here is completely general and can be applied to 
many other physical systems. 
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FIG. 2: Evolution of the state along the no-jump trajectory on 
the Bloch sphere, for the spontaneous decay case (ri = \a z , 
T2 = q<t_). The evolution is a smooth spiral converging to 
the lower state: while the state rotates about the z axis, it 
is smoothly brought towards the lower state with a velocity 
given by the spontaneous decay rate a. The geometric phase 
is given by the area enclosed in the path shown, where the 
last segment is the geodesic connecting initial and final point 
of the evolution. 
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